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Abstract 

We evaluate the axial anomaly for the general Ginsparg- Wilson fermion operator 
D = -Dc(I + RDc)~^ with R = rl. For any chirally symmetric Dc which in the 
free fermion limit is free of species doubling and behaves like i'y^p^ as p ^ 0, the 
axial anomaly tr[75(i?L')(x, x)] for U{1) lattice gauge theory with single fermion 

2 

flavor is equal to ■^^e^j,,^XfjF^jy{x)F\^{x + jl + v) up to higher order terms and/or 
non-perturbative contributions. The FF term is r-invariant and has the correct 
continuum limit. 
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1 Introduction 



In 1981, Ginsparg and Wilson [1] formulated a criterion for breaking the chiral 
symmetry of the massless Dirac operator on the lattice. 



D75 + 75^ = 2aD^^RD (1) 

where R is any invertible hermitian operator which is local in the position 
space and trivial in the Dirac space, and a is the lattice spacing. The underlying 
reason for breaking the continuum chiral symmetry on the lattice is due to 
the Nielson-Ninomiya no-go theorem [2] which states that any Dirac operator 
on the lattice cannot simultaneously possess locality, free of species doubling 
and the chiral symmetry. However, one prefers to violate the chiral symmetry 
rather than the other two basic properties, is due to the fact that if the chiral 
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symmetry breaking is specified by the RHS of (1) having one 75 sandwiched by 
two Dirac operators, then not only the usual chiral symmetry can be recovered 
in the continuum limit a — > 0, but (1) also gaurantccs the remnant chiral 
symmetry on the lattice, that is, the action A — ipDip is invariant under the 
finite chiral transformation on the lattice 

ijj exp[075(I - aRD)]'ijj (2) 
^ -0 exp[^(I - aDR)^5\ (3) 

where ^ is a global parameter. The infinitesimal form of (2) and (3) was first 
observed by Liischer [3]. In fact, the GW relation (1) can be generahzed to 
accommodate the asymmetric finite chiral transformations on the lattice [4] 
having two different hermitian operators, say S and T in (2) and (3) respec- 
tively, by replacing 2R in (1) by S + T. Furthermore, the particular form 
of chiral symmetry breaking on the RHS of (1) also implies the chiral Ward 
identities, non-renormalization of vector and flavor non-singlet axial vector 
currents, and non- mixing of operators in different chiral representations [5] . 

In general, on the lattice, given any chirally symmetric Dirac operator Dc 

^c75 + 75^c = 0, (4) 

which is free of species doubhng but nonlocal as a consequence of the Nielson- 
Ninomiya no-go theorem, the general solution [6] 



D ^ D^(l + aRD^)-^ (5) 

is a chiral symmetry breaking transformation which gives a class of Dirac op- 
erators satisfying the GW relation (1). The general solution (5) also implies 
that any zero mode of D is also a zero mode of D^., and vice versa. That is, the 
zero modes of D are i?-invariant. Then the index of D is also /^-invariant and 
is equal to the index of Dc- Thus the chiral symmetry breaking transformation 
(5) cannot generate a non-zero index for D if the index of is zero [7]. It 
has become clear that a basic attribute of D{D^ should be introduced, and 
it is called topological characteristics in ref. [7,8]. In general, the topological 
characteristics of a Dirac operator cannot be revealed by any perturbation 
calculations. Therefore, if we obtain non-zero axial anomaly for D in a per- 
turbation calculation, it does not necessarily imply that the index of D must 
be non-zero for topologically non-trivial background gauge fields. The most 
reliable way to determine the topological characteristics of a Dirac operator 
is to perform the following numerical test [9,8]: turn on a topologically non- 
trivial smooth background gauge field with constant field tensors, then check 
whether there are any zero modes of D, and measure the index of D versus 
the topological charge of the background gauge field. It should be emphasized 
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that the topological characteristics of D is a basic attribute of D, which is 
not due to finite size effects, but persists for any lattice sizes, at any lattice 
spacings and in any background gauge configurations. 

In Ginsparg and Wilson's original paper [1], the axial anomaly for any D 
satisfying Eq. (1) was derived, and it agrees with the continuum axial anomaly 
if D is free of species doubling and in the free fermion limit behaves like 
^ P ~^ 0- However, in their derivation, Eq. (1) is used to eliminate 
R in some of the intermediate expressions, then at the final steps of their 
computations, D is replaced by Dc, which is equivalent to setting R = 0. 
Strictly speaking, their procedures are not completely self-consistent since R 
should be kept nonzero throughout the entire computation, otherwise Eq. 
(1) cannot be used to eliminate R from any expressions containing D'^^RD. 
In other words, setting i? = at their final steps of integrations actually 
invalidates their previous steps of using Eq. (1) to eliminate R. The proper 
procedure should keep i? 7^ throughout the entire calculation, and then 
shows that the axial anomaly is independent of i?, finally the limit i? = can 
be safely taken. This motivated us to re-derive the axial anomaly for a general 
D satisfying the Ginsparg- Wilson relation, with the proper procedure, and in 
the context of recent developments. Furthermore, the realization of the Atiyah- 
Singer index theorem on a finite ( infinite ) lattice relies on the perturbation 
result of the axial anomaly as well as the topological characteristics of D. This 
provided us further impetus to carry out the tedious computations and present 
the details of our derivation in this paper. We note in passing that in ref. [10,11] 
the chiral anomaly for the overlap-Dirac operator [12,13] is calculated, and 
their results agree with the continuum axial anomaly. However, the overlap- 
Dirac operator is only one of the solutions satisfying the Ginsparg- Wilson 
relation, and of course their results do not imply that the same axial anomaly 
will be obtained for a general D with another Dc. 

If we also require D to satisfy the hermiticity condition 



L>t = ^^D-f^ (6) 

then Dc also satisfies this condition, and becomes an antihermitian operator 
{ D\ — —Dc ) which has one to one correspondence to a unitrary operator V , 

D. = ^ (7) 



where V also satisfies the hermiticity condition 75F75 = . Thus the inverse 
operator in Eq. (5) must exist, and the general solution D is well defined. 

In this paper, we evaluate the axial anomaly ii[y^{RD){x,x)] for R diagonal 
in the position space, i.e., R — rl with parameter r. Then from Eq. (5), we 
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have 



D^D^{l + rD,) 



-1 



(8) 



We shall restrict our discussions to the U{1) gauge theory with single fermion 
flavor. However, it is straightforward to generalize our derivations to lattice 
QCD. For any D^. which in free fermion limit is free of species doubling and 
behaves like i'Ji^Pn as p — > 0, our perturbation calculation shows that the FF 
term of the axial anomaly tT[y5{RD){x,x)] is independent of r and has the 
correct continuum limit, i.e.. 



r tr[75i:>(a;, x)] = ■^^€,^^xaF^,uix)Fxaix + fi + u)+ other terms (9) 
where the field tensor on the lattice is defined as 



The other terms in (9) in principle cannot be computed directly by any pertur- 
bation calculation to a finite order, however, their sum over the entire lattice 
can be determined and has significant impacts to the index theorem on the 
lattice, as we will show in section 2. 

The outhne of this paper is as follows. In section 2, we derive the divergence of 
the axial vector current for the Dirac operator satisfying the general Ginsparg- 
Wilson relation, and to set up the theoretical framework for the perturbation 
calculation in section 3. The topological characteristics of D is shown to emerge 
naturally as an integer functional of D, after the axial anomaly is summed over 
the entire lattice. In section 3, we perform the derivation of the axial anomaly. 
In section 4, we conclude and discuss. In appendix A, we derive an identity for 
the FF terms. In appendix B, we derive some useful properties of the kernel 
for the vector current which are used in the derivation of axial anomaly in 
section 3. In appendix C, we prove an identity for the Ginsparg- Wilson kernel 
of the vector current. 



2 The cixial vector current and its divergence 

In this section we derive the divergence of the axial vector current A, D) 

for the Dirac operator D satisfying the Ginsparg- Wilson relation of the general 
form [4] 




(10) 
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L>75(I - SD) + (I - DT)-f5D = (11) 

where 5" and T are arbitrary invertible hermitian operators which are local 
in the position space and trivial in the Dirac space. The action for exactly 
massless fermion in a background gauge field is 



A^J2'^M^,y;A)il;y (12) 

x,y 

where X and y are site indices, and the Dirac indices are suppressed. Then the 
action A is invariant under the chiral transformation 



^ ^ cxp[d-,'5(l- SD)]tP (13) 
^ V exp[^(I - L>r)75] (14) 

where ^ is a global parameter. Hence, the divergence of the ( associated 
Noether current ) axial vector current, d^J^{x), can be extracted from the 
change of the action 5 A under the infinitesimal local chiral transformation at 
the site x, 

^^. + 5^.75[(lI-5/^)^^]x (15) 
^^^^^ + 5eSiJ-DT)],l, (16) 

with the prescription 

A^ A + Se^di'Jlix) (17) 
Then we obtain 

-{iPDT),j,{D^), - {i,D),^^{SDxl^), (18) 
which satisfies the conservation law 

Y.d'Jl{x)^Q (19) 

X 

due to the exact chiral symmetry (11) on the lattice. Now we take the lattice 
to be finite and with periodic boundary conditions, then we define 5^J^(a;) by 
the backward difference of the axial vector current 
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such that it is parity even under the parity transformation, and the conserva- 
tion law Eq. (19) is also satisfied. 



To evaluate the fermionic average of the divergence of the axial vector current 
in a fixed background gauge field, 



{di}\{di}\d^ Jl{x) exp{-iljDilj) 
J [dV][#]exp(-'0DV) 



(21) 
(22) 



one would encounter which is not well defined for the exactly massless 
fermion in topologically nontrivial gauge background. Thus, one needs to in- 
troduce an infinitesimal mass m which couples to the chirally symmetric Dirac 
operator Dc in the following way [4] 



D = {Dc + m) 



1+-{S + T)D, 



(23) 



and then evaluate the fermionic average (21) with D replaced by D, and finally 
take the limit ( m — > ), i.e.. 



■ lim — 



[dij][dilj]d''jf^{x) exp{-iljDilj) 
Z^J [#][#] ex^{-i^Di)) 
Substituting d^J^{x) by Eq. (18) and using Eq. (23), we obtain 



(24) 
(25) 



d'^r^ix)) = -tr[(I - L'T),,75] - tr[(I - SD),,^,] 



-|- m tr 



1--{S + T)D 



D-\1-DT)] 75 

J XX 



+ m tr N (I - SD)D 



-1 



1--(S + T)D 



75 



(26) 



where tr denotes the trace in the Dirac space. The first two terms on the RHS 
of (26) is equal to 



tr{^,[{S + T)DU} 
while the last two terms in the limit ( m — > ) give 



(27) 



2 m trfcSs] 



(28) 
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which can be rewritten as 



^^^^tm^ML (29) 

a,/3 s m + Xs 

where (ps and As are normahzed eigenf unction and eigenvalue of D, 



y (3 

EE[C(^)]>:'(^) = <^..' (30) 

X a 

In the hmit ( m — > ), only zero modes of D contribute to (29) and the result 
is 

-2YM^;ix)]^cp;{x) (31) 

s=l t=l 

where and (p^ are normalized eigenfunctions oi D { ) with eigenvalues 
Xg — Xt — and chiralities +1 and —1 respectively. Therefore in the limit 
m — > 0, Eq. (26) becomes 



d^4{x))^tr{j,[(S + T)DU} 

+2j:[ci>t{x)]Uti^) - 2E[0r(^)]Vr(^) (32) 



s=l t=l 



This is the anomaly equation for D satisfying the general Ginsparg- Wilson 
relation (11) which is the exact chiral symmetry on the lattice, where the 
axial vector current J^{x) is the associated Noether current. 

On the other hand, as usual, if one considers the action built from the chirally 
symmetric part of D [1] , 



A = ^(^-75^^75) V' (33) 

then As has the usual chiral symmetry, and the divergence of the associated 
Noether current can be extracted from the change of the action 6 As under the 
infinitesimal local chiral transformation at the site x 
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with the prescription (17), and one obtains 

~[iPD{S + T)^5DU, - ^iPx[D^5{S + T)D'iPl (34) 

which also satisfies the conservation of total chiral charge, Eq. (19), due to the 
usual chiral symmetry. Although Eq. (34) looks fairly different from Eq. (18), 
however, the fermionic average of (34) in a background gauge field [14] is equal 
to that of (18), i.e., Eq. (32). It is evident that the difference between (18) and 
(34) has no physical consequences. In fact, it is possible to redefine d^J^{x) 
of Eq. (18) in many different ways provided that it satisfies the conservation 
law, Eq. (19), and its fermionic average agrees with Eq. (32), however, the 
corresponding J^{x) is not equal to the Noether current associated to the 
exact chiral symmetry on the lattice (11). For example, if one redefines (18) 
as 

a'^jj(x) =^/;,75(D^). + (^I^).75^. 

-'4)^{DT-'ir,Dij)^ - {i)D'-fr^SD)^i)^ 
= [^^75(1 - SD%^^ + V.[(I - DT)j5Dijl 
= [iPDj,{l - SD)l,p, - iP,[Dj,{l - SD),pl (35) 

where Eq. (11) has been used in the last equality, then the fermionic average 
of Eq. (35) obviously agrees with Eq. (32), and the conservation law, Eq. (19) 
is also satisfied. 

Likewise, the divergence of the vector current can be extracted from the change 
of the action 5 A under the infinitesimal local transformation 

with the prescription 
then we obtain 

dnJn{x) = {i)D)^'^^ - '4)^{Dil))^ (36) 
where d^J^^x) is defined by the backward difference of the vector current 
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(37) 



such that it is parity even under the parity transformation, and the conserva- 
tion law due to the Uv{l) symmetry, 



J2dMx)^0 (38) 

X 

is satisfied on a finite lattice with periodic boundary conditions. If the vector 
current is expressed in terms of the kernel Kf^ as 

Jf,{x; A,D)^Y, i^x+yK^ix, y, z\ A, D)ip^+^ (39) 

then by comparing Eqs. (34) and (36), the axial vector current satisfying (34) 
can be written as 



4(x; AD)^Y1 i^-+yKl{x, y, z; A, L>)V'.+. (40) 
y,z 

where the kernel K^^ is related to by 



Kl^\{K,^,-^,K,) (41) 

We note in passing that for the JT"^ defined in (35), one could not find a simple 
relationship between the kernel of this axial vector current and the kernal 
of the vector current, since (35) is non-linear in D. 

Now summing Eq. (32) over all sites of the lattice, the LHS is zero due to the 
conservation law (19), then the RHS gives the so called index theorem on the 
lattice [14,3,12] 



7V_ - iV+ = ^ {75[(^ + T)DU} (42) 

^ X 

where N^{N_) denotes the number of zero modes of positive ( negative ) 
chirality. It has been shown in ref. [7] that the index is invariant with respect 
to S and T and is equal to that at the chiral hmit. 



7V_ - iV+ = ^hm^ i E tr {^,[{S + T)DU,} ^ ^Yl t^N^.,.] (43) 
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where Eq. (7) has been used. We note that a priori, there is no compeUing 
reasons to gaurantee that D has zero modes in topologically non-trivial sectors. 
It could happen that D is local and free of species doubling in the free fermion 
limit, but turns out to have zero index in any background gauge fields [8]. 
In that case, D is called topologically trivial, and the index theorem (42) is 
trivially satisfied with both sides equal to zero, however, it does not correspond 
to the Atiyah-Singer index theorem in continuum. Presumably the index of 
D is a topological and non-perturbative quantity, therefore the topological 
characteristics of D cannot be revealed by any perturbation calculculations 
at finite orders. We refer to ref. [7,8] for further discussions on topological 
characteristics of D. As we will show later in this section, the topological 
characteristics of D emerges naturally as an integer functional of D, after the 
axial anomaly is summed over all lattice sites. 

In the next section, we will show that the first term on the RHS of Eq. (32), 

tr{75[(5 + T)D],,,} (44) 

reproduces the topological charge density 

p{x) = Y^ef,yXaFi,^{x)Fxa{x + fi + i)) (45) 

up to higher order terms and/or non-perturbative contributions. The field 
tensor F^,y{x) on the lattice is defined in (10) for QED. We note that in (45) 
the positions of the field tensors F^,/ and Fxa are chosen at x and x + fi + u 
respectively such that p{x) satisfies 

EMa^) = (46) 



for any local deformations of the gauge field [15]. In the continuum limit, p{x) 
agrees with the Chern-Pontryagin density in continuum. 

In order to extract the term which is quadratic in gauge fields from (44), we 
consider the following operator H [14] which is defined as 



- E-.I/atT^tA:, \0{x)] (47) 



'SAM SA^iy) 



where 0{x) is any observable. It is evident that the operator Ti. picks up the 
terms quadratic in gauge fields from the observable O and converts them into 
a constant. 

The operator H acting on (45) gives 
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n[p\ 



167r2 



^l3-isFaf3{x)F^s{x + a + P) 



A=0 



(4J 



(49) 



where Eq. (10) has been used. We note that it is not necessary to take the 
hmit Afj_ = after the differentiations since p is quadratic in gauge fields. The 
proof of Eq. (49) is given in appendix A. We note in passing that 

where the operator Ti' is similar to Ti. defined in (47) but without imposing 
Aj^ = after differentiations with respect to the gauge fields. However we 
exclude the case of p having both field tensors located at x since it does not 
satisfy Eq. (46). 

On the other hand, if we have Ti act on (44) and obtain 



n{tr{j5[iS + T)Dl,})^ 



(50) 



then we can infer that 



tr{75[(5 + T)D]^^} = —-^e^,xaF^,u{x)Fxa{x + p + u) + other terms (51) 

where "other terms" denotes those terms which cannot be determined by 
the second order perturbation calculation using Ti, which in general con- 
sists of higher order terms in A^ and/or derivatives, plus terms due to non- 
perturbative and/or topologcial effects (if any). The other terms in (51) in 
general cannot be computed directly by any perturbation calculation to a fi- 
nite order, however, their sum over the entire lattice can be determined and 
might have significant impacts to the index theorem on the lattice, as we will 
show later. 

Although Eq. (51) refers to the infinite lattice, however, we expect that it 
also holds for finite lattices. The argument [7] is as follows. If D is local, the 
boundary effects enter as ~ exp{—m{R)L/a), where L is the lattice size, a is 
the lattice spacing, m{R) is a monotonic increasing function of i? = rl ( take 
the simplest case ) with m(0) equal to zero ( this is equivalent to that Dc is 
nonlocal ) and m(oo) a positive constant. As L — oo, the finite size effects 
vanish and the axial anomaly is given by Eq. (51). For L is finite, one might 
naively expect that the L dependence would enter the FF term and "other 
terms", through the variable m{R)L/a. However, the FF term cannot have 
it! dependence, otherwise it would be in contrary to the fact that the index 
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of D is i?-invariant [7]. Consequently, if L is gradually decreased from infinity 
toward a finite value, all dependence of L and R only resides in "other terms" . 
So, Eq. (51) also holds for finite lattices provided that D is local. 

After summing Eq. (51) over all sites of the lattice, we have 
2{N_ -N+)=Y: tr {75[(^ + T)Dl,,} 

X 

= H W(^-^M'^(^)-^Aa(2; + /i + i>) + Xl( Other terms ) (52) 



Now we consider the topologically nontrivial background U{1) gauge field on 
a 4-dimensional torus ( e [0, L^j], /i — 1, ■ ■ ■ ,A ), 

eAW = H^-H^ + 4»).in(H^..) (53) 
e^,M = ?0 + 4")sm(H^x,) (54) 



eM.) = ?^-^Sp + ArsiJ?^.A (55) 

e^4(x) = ^ + 4°^sin(^X3) (56) 



where qi, q2, rii, - • ■ ,n4 are integers. The global part is characterized by the 
topological charge 



-.2 



Q = j d^x e^^xa Fi^u{x)Fxa{x) = qiq2 (57) 



which must be an integer. The harmonic parts are parameterized by four 

real constants hi. h2, /13 and /i4. The local parts are chosen to be sinusoidal 
fiuctuations with amplitudes Ai'\ ^2°"*) ^3 ^^'^ A^^\ and frequencies ^f^, 
^ and ^ respectively. The discontinuity of Ai(a:) ( A^ix) ) at 2:2 = L2 
( 0:4 = I/4 ) due to the global part only amounts to a gauge transformation. 
The field tensors F12 and F34 are continuous on the torus, while other F's 
are zero. To transcribe the continuum gauge fields to link variables on a finite 
lattice with periodic boundary conditions, we take the lattice sites at x^ — 
0, a, (A'^^ — l)a, where a is the lattice spacing and — N^a is the lattice 
size. Then the link variables are 



Ui{x) = exp [ie^i(a;)a] (58) 

r 27iqiXi' 
L'2(x)=exp ieA2{x)a + idx^^(^N2-i)a- 



(59) 
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Usix) = exp [ieA3(x)a] (60) 
U4{x) =exp 



ieA4(a;)a + i(5^4,(Ar4_i)a^^Y^ (61) 



The last term in the exponent of U2{x) ( U4{x) ) is inchidcd to ensure that 
the field tensor F12 ( -F34 ) which is defined by the ordered product of hnk 
variables around a plaquette [ Eq. (A. 8) ] is continuous on the torus. Then the 
topological charge of this gauge configuration on the finite lattice is 



Q = Yl ^fii^X'TF^,„{x)Fxa{x + jl + i>)^ 51^2 (62) 

which agrees with the topological charge [ Eq. (57) ] on the 4-dimensional 
torus. From Eq. (52), we obtain 



^( other terms ) = 2(A^_ - N+ - giga) (63) 

X 

which is also an integer. Since D does not have any genuine zeromodes in the 
topologically trivial gauge background, thus the integer (63) must be propor- 
tional to qiq2 and can be represented as 

other terms ) 

X 

= {c[D] -l)Y, ■Y^(^ni'\aF^Ax)Fxa{x + fl + D) 
where c[D] is an integer functional of D. Then Eq. (52) becomes 



2 

^^e^^XaF^,^.{x)Fxa{x + /f + i>) 

^c[D] Q (64) 

where c[D] = c[Dc\ due to the invariance of the zero modes and the index under 
the chiral symmetry breaking transformation (5) [7]. In general, Eq. (64) holds 
for any smooth background gauge configurations on a finite ( infinite ) lattice 
provided that the topological charge on the lattice (62) is an integer and the 
axial anomaly on the lattice satisfies Eq. (50). It is remarkable that, if any of 
these higher order, non-perturbative and topological contributions to the axial 
anomaly does exist, then their total effects to the index can only enter as an 
integer ( c[D] — 1 ) multiple of the topological charge of the background gauge 
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field. Wc emphasize that the emergence of c[D] is not due to the finite size of 
the lattice but an intrinsic characteristics of in particular, when c[D] = 0, 
then the axial anomaly vanishes at each site of the lattice, independent of the 
size of the lattice. In general, the topological characteristics, c[D], is an integer 
functional of D, which in turn depends on some parameters of D as well as the 
gauge configuration. c[D] could become chaotic when the background gauge 
field is rough, as first demonstrated in ref. [8] . For smooth gauge configurations, 
D can be classified according to its topological characteristics as follows [7]. 
If c[D] = 1, then D is called topologically proper, else if c[D] = 0, then D 
is called topologically trivial; else c[D] — integer ^ 0, 1, then D is called 
topologically improper. Only for D is topologically proper, i.e., c{D\ = 1, Eq. 
(64) can realize the Atiyah-Singer index theorem on a finite ( infinite ) lattice. 
[ Note the invariance of the topological charge in Eqs. (62) and (57) ]. Eq. 
(64) provides the theoretical understanding of the numerical results [9,8] for 
exact zero modes satisfying the Atiyah-Singer index theorem even on very 
small lattices in two dimensions as well as in four dimensions. 

We note that in ref. [8], using the exact refiection symmetry and the exact 
solution of the free fermion propagator, the ( lowest order ) perturbation 
theory is shown to break down at the topological phase boundaries in the tuq 
parameter space of the overlap-Dirac operator. The zero index of D at the 
phase boundaries can be interpreted as c[D] = 0, due to non-perturbative 
and/or topological contributions. In general, c[D] incorporates all kinds of 
contributions from all fermionic modes. 

We also note that in the operator Ti., the gauge fields are set to zero after the 
differentiations with respect to the gauge fields. This implies that only free 
fermion propagators are needed in the evaluation of the LHS of Eq. (50). In 
the next section, we will show that Eq. (50) is indeed satisfied by the Ginsparg- 
Wilson Dirac operator (5) with R — rl, provided that Dc in the free fermion 
limit is free of species doubling and behaves like i'JuP^ as p — > 0. 

Recently Liischer [15] proved that for t/(l) lattice gauge theory, and for S+T — 
2, if the axial anomaly q{x; A, D) = trlj^Dx^x] satisfies 



(65) 



X 



for any local deformations of the gauge field, then 



q{x; A,D) = j eij,„xaF,j,u{x)Fxa{,x + fi + u) + d^G^{x; A, D) 



(66) 



where 7 is a constant times an integer factor, and 
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d^G^ix; A,D)^Y1 [Gf^ix; A, D) - G^{x - fi; A, D)] . (67) 

The explicit form of the current Gfj,{x; A, D) is supposed to be very comph- 
cated. As discussed in ref. [7], Eq. (66) can be generahzed to any D satisfying 
(11) with all the S and T dependences residing in the term d^Gfj,{x; A, D), 
while the FF term depends on the topological characteristics c[D] which is 
(,5, r)-invariant ( c[D] = c[Dc] ), i-e., 



q{x; A, D) = \iT{^^[{S + r)D],,,} (68) 
= V c[D^]e^,xcF^,{x)Fxa{x + /t + i>) + d^G^{x- A, D) (69) 

where the current G^^x] A, D) in general is a functional of S and T. Although 
Eq. (66) is proved for the infinite lattice in ref. [15], however, as discussed in 
ref. [7], if D is local, then the boundary effects which depend on {S,T) and 
the lattice size L, can only enter the term d^Gf^{x; A, D), thus Eqs. (66) and 

(69) are also true for finite lattices. The same argument has been presented 
in details to assert that Eq. (51) also holds for finite lattices. [ See the second 
paragraph after Eq. (51). ] 

Now applying the operator 7i to Eq. (69) and using Eq. (49), we obtain 

H [q{x; A, D)] ^-f' H c[D] eapjsFap{x)F^s{x + a + $) 
+n[dA{x;A,D)] 
= 87' e^^Aa + n[d^Ga(x; A, D)] (70) 

where c[D] is presumably non-perturbativc and thus cannot be determined 
by the second order operator TC acting on q{x;A,D). Hence D can only be 
assumed to be topologically proper and c[D] is set to unity throughout the 
entire perturbation calculation. Then the assertion of Eq. (50) implies that 
the LHS of Eq. (70) is (S*, T)-invariant. On the other hand, on the RHS of 

(70) , the first term is {S, T)-invariant, but the second term in general depends 
on S and T. This implies that H[daGa{x; A, D)] must vanish identically, and 
Eq. (70) gives 



and (69) becomes 



q{x; A, D) = ^ c[D^]e^,xaF^,,{x)Fxa{x + /i + i>) + d^,G^{x- A, D) (72) 
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For the U{1) gauge fields defined in Eqs. (58)-(61), summing Eq. (72) over all 
lattice sites yields 

N_-N+ = c[D] Q 

in agreement with Eq. (64). 

It is interesting to note that the integer factor in the 7 of Eq. (66) turns out can 
be used to accomodate non-perturbative and/or topological effects, and it is 
identified to be the topological characteristics, c[D] which is first introduced in 
ref. [7,8] in the study of the index of D with rcpscct to the background gauge 
fields. It also provides plausible explanations to some seemingly paradoxial 
situations that one obtains the correct axial anomaly in the perturbation cal- 
culation for a given D but the numerical ( nonperturbative ) computations 
give exactly zero axial anomaly at each site, for any lattice sizes and for any 
gauge configurations. 



3 The Axial Anomaly 

In this section, we assert Eq. (50) for 5" + T = 2rl by evaluating 

l2 = n[2 r tr(75D.,.)] 

where D satisfies the Ginsparg- Wilson relation (1) with it! = rl, i.e.. 



(73) 



D^D,(l + rD,)-^ (74) 

for any chirally symmetric Dirac operator Dc which in the free fermion limit, 
is free of species doubling and behaves like i^y^p^ as p ^ 0. 

In general, we can write Dc in the momentum space as 



Dcip) = il,C,{p) (75) 

where C^{p) are arbitrary functions which satisfy the following properties in 
the free fermion limit : 

(i) C^,{p) ^ p^ as p ^ 0. 

(ii) C^{p) has no zeros in the entire Brillouin zone except at the origin p = 0. 
Then using Eqs. (74) and (75), we obtain 
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^ 1 -iP 



(76) 
(77) 



First we perform the differentiations with respect to the gauge fields. The 
general formula is 



S 5 



{TiU)) 
1 



5A,{Q) 6A^{m) Z(U) 



T{U)e 



-'il>D{U)'il) 



12 



(78) 



i=l 



where U denotes link variables and T{U) is an arbitrary functional of link 
variables and fermion fields, and 



P2 
P3 

Pa 
P5 
Pe 
P7 
Ps 
P9 



SA^{m) 



6 



\^SA,{0)6A^{m 

' :TiU)]U ' 



6A,{oy 

-D{U)ij 



SAf,{m) 



D{uyijj 



5A4oy 



\\5A^{m) 
' S S 



5A,(oy 



T{U) 



^SAij,{m 



y5Ai,{m) 
-TiU) 



5A,{Q) 



D{U)i) 



5 



M,(0) 



(^(^)(^M?(o)'^H>< 



D{U)i) 



y, 



dA^{m 



-D{U)7jj 



D{U)ip 



5A^{m) 



SAM 



11 



{T{uy y 



5A,{Q) 5A^{m 



-D{U)y 



79 
80 
81 
82 
83 
84 
85 
86 
87 



89 
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To evaluate I2 in (73), we set 

T{U) = {tljDR^,D)niJn = r{^ljD^5D)niJn (91) 

Then (T(f/)) = — r tr[75D„„] which vanishes in the free fermion hmit. So, 
Pio, Pn and P12 do not contribute to l2- Next we express the vector current 
in terms of the kernel as 

However, the vector current satisfying the divergence condition Eq. (36) is 
not unique. A general and explicit realization is given by Ginsparg and Wilson 
[1] with K^{k, U) equal to D^+i^k+jiU) times the sign of {i— j)^ and times 
the fraction of the shortest length paths from k+ j to k + i which pass through 
the hnk from k to k+fi. The Ginsparg- Wilson kernel can be expressed in terms 
of the derivative of the Dirac operator with respect to the gauge field. 



n^"'SA,{k) 



Dmn{U)lPn (92) 



where the second equality is proved in appendix C. In the free fermion limit, 
the action is translational invariant, — D^_n, then Eq. (92) becomes 



•^m(^) = Z] ipk+iK^{i, i - j)'il^k+i-j = i X] V'n 



SA^ik) 



Dmn{U)lpn 



(93) 



A=0 



with 



i-j)^ sign(j^) f^{i, j) Dj (94) 

where fn{i,j) is equal to the fraction of the shortest length paths from to j 
which pass through the link from i — (ito i. 

We note that Hasenfratz [16] also showed that 

m,n ^^nKf^) 

satisfies the divergence condition, Eq. (36). However, an explicit realization of 
the kernel has not been prescribed. Furthermore, in the free fermion limit, the 
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kernel is shown [16] to satisfy the so cahed sum rules which are equivalent to 
Eqs. (B.8) and (B.9) proved in appendix B as well as in ref. [1]. 



According to the properties of K^i^k, i,j;U) discussed above and the following 
simple identity of the derivative of a link variable with respect to the gauge 
field, 



SA^{n) 



(95) 



we obtain 



K,{k,z,j;U) 



Cnu{n)Ki^{i,i - j) 



(96) 



A=0 



where c^i,{n) are some constants independent of the gauge fields. Some useful 
properties of Ki^{i,i — j) are derived in Appendix B and will be used in our 
evaluation of l2- 

Using Eqs. (91), (92), (93), and (96), the expressions of Pis in Eqs.(79)-(87) 
become 



Pi^-{{'ij^DR-f5. 
P2 = ici,u{0) {{^DR'-ir,D)ni)nJii{m 



P3 = t({ j^j-^ii^DR^.D)^,^;, 



'n 



P4^i 



----Ai^DR^^D)r,i,, 
6A^{m) 

S 



A=0/ 



'n 



Pe = i^{^DR-i^D\i)^J^{m)) ( J,(0)) 
\dA^{m) 



MO)] 



A=o/ 



A=Oi 



(97) 
(98) 

(99) 

(100) 

(101) 
(102) 
(103) 
(104) 
(105) 



n = -i (| ^^^(^L>P75^)nV'n ^ (Urn)) 
where the free fermion limit has been imposed for the RHS of Eq. (73), i.e. 
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9 

/2 = -2 E{E^-^APi} (106) 

i=l m,n 

In above expressions, only Pi has non- vanishing contributions to l2- This can 
be shown in the following. First consider the common factor (Ji/(0)) of Pq and 
P7. Using Eqs. (93) and (94), we obtain 

= ^tr[D:]sign(j,) /(z,-j) D_j] 



-tr(l)^sign(j,)/(^,j)=0 (107) 



since J2if{'i,j) = J2if{h —j)- Similarly, we can show that {J^{m)) = since 
it is translational invariant. Then Pg, P7, Ps and Pg in Eqs. (102)-(105) are 
zero. 

For P2, the factor ({;ipDR'^^D)nipnJn{'<^)/ enters (106) to give 



namx {{'^DR'yr,D)n'^nJn{m 

m,n 
m,n,i,j 

m,n,i,j 

-tr[^5Do] tr[K^{i,i-j)Dz]]} 

= 

where the identity 



tr(75) = tr(757^) = tr(757^7^) = tr(757^7^7^) = (108) 

has been used, and the fact that each factor of D, D^^ or K^j can give at most 
one gamma matrix in the free fermion limit [ Eqs. (76)-(77) and (94) ]. 

For P5, it enters (106) to give 



1 / 5 S 



A=0/ 

' ) (109) 



A=0/ 
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where Eq. (1) has been used. Using Eq. (95), we immediately see that this 
expression only involves trace operations on terms containing one 75 and less 
than four 7^5 matrices, thus it must be zero. 

For P3, it enters (106) to give 



Jn{m) 



A=0/ 



2 / _ 

m,n,k,i,j 



Tn,n,k,i,j 



[D„Xi_j_k[Ku{k, k - n)75 + -i^K^{k, k - n)]D,^\_^_^K^{i, i - j 

-tr {Dzl,[K{k, k - n)75 + ^sKik, k - n)]) tr {DzjK^{i, i - j)) } 
= 

where Eqs. (1) and (93) have been used in the first and the second equalities. 
From Eq. (94), we see that [K^^^k, k — n)^^ +'~^^Ky{k, k — n)] is proportional to 
75, then using the trace identity Eq. (108) to give zero in the above expression. 
Similarly, we show that P4 also enters (106) to give zero. 

Finally only Pi remains in Eq. (106), i.e.. 



/2 = -2 ^ n^mxPi 

m,n 



= Yl n„mx{:ipk{Dk-nlb + '^-,Dk-n)lpn ■ 
m,n,k,i,j,s,t 

'ipm+iK^ii, i - j)i/jm+i-jAK^{s, s - t)i^s-t) (110) 

where Eqs. (1) and (93) have been used. After the fermion fields are contracted, 
Eq. (110) becomes 

h^Ia + h + h + h (111) 

where 
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4 = - E n,rnMl5D-l^_,K,ii,i-j)D-\,_._^K,is,n)] (112) 



m,n,i,jyS 



m,n,i,s,t 

D-Xi-,-sKA^,^-t)D7lt-k] (114) 

-^d = - E "'fT"^Atr[75-Dfc-n-D~^5Xi.(s, s-t) 
m,n,k,i,j,s,t 

D:1,_^_^K,{z, I - j)D-\^_^_,] (115) 

where those terms involving trace operation on products of one 75 and less 
than four 7^5 matrices are zero and have been dropped. Using the Fourier 
transforms 



■—7T 

K,{m,m-l) = J 0^^e^P-'e-^^'<-^-^^K,{p,p + p') (117) 

— TV — TT 

and the identity 



(118) 

we obtain 



Jp4im+j-s) ip^s-ipen 

m,n,l,j,s ^p^y 

tr[^5D-\pi)K^{p,,p,)D-\p,)K,ip,,pe)] 

= J 5{Pa)8{Pb)dx^d^a^ j tr[75^"'(p + Pb)K^{p + Pb:P + Pa+ Pb) 

PaPb V 

D-\p + Pa+Pb)K,ip + Pa+Pb,p)] (119) 

The differentiation with respect to Pa and pb will generate many terms in 
the last expression. However, due to the identity (108), only those terms con- 
taining the product O571O27374 and its permutations can have nonzero con- 
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tributions, thus the final result is proportional to tr(757^7,^7A7cr) = 4:eij,i,xa- 
Then it is obvious to see that those terms containing d^^^ d'^^ D~^{p + Pa + 
Pb), or d^^^ dj^'' Ki,{p + Pa + Pb,p) must vanish since they are symmetric in A 
and a. Furthermore, the terms containing d^^^K^, d'^^K^^ d^^^ K^^ dj^^'Ki,, or 

d^^^ d^^ Kn{p+pi„p+Pa+Pb) would become zero after integrations over pa and 
Pb with the delta functions in (119), due to the following identities ( proved 
in Appendix B ) 



d',K^{p,P + P%'=o^(^ld^d, - ^S^,d,^ D{p) (120) 
d',K,ip + p',p)\p,=o=[ld,d, + ^d,xd,^ Dip) (121) 

which are symmetric in /j, and A. Hence, the only nonzero term in (119) is 

4 = / tr{j5[daD-\p)]K,{p,p)[d,D-'{p)]K,{p,p)) 

\r{^,[d,D-'{p)][d,D{p)][dxD-'{p)][d,D{p)]) (122) 



where the identity ( proved in Appendix B ) 

K,{p,p)=td,D{p) (123) 

has been used. By performing the same analysis on If. and la, we obtain 

m,n,l,s,t 



tT[^,D-\pi)K,ip2,P3)D-'ip^)K^ip5,pe)] 
tT[^,D-\pi)K,{p,,Ps)D-\ps)K^{ps,Pe)] 

S{Pa)S{pb)di^^d!^^ I tr[75D-^(p +Pa+ Pb)K,{p + Pa+ Pb:P) 



PaPb 



D-\p)K^{p,p+Pa)] 

= (124) 

where the zero is essentially due to the presence of D~^{p) which does not 
depend on pa or ph. 



23 



m,n,k,l,j,s,t ^p^y 

e'P'^'-'hr[^5D{pi)D-\p,)K,{p,,p,)D~\p,)K,{pe^ 

= - ^ n^mx J er^V^^e^Vi{n-va)^p^ra 

Pl,P2,P4 

tr[75L'(pi)L'-^(p2)i^M(P2,P4)£'-'(p4)i^.(p4,Pl)^-'(pi)] 

= / 5(p„)5(p,)9f9W|tr[75L'(p)L'-^(p + p,) 

PaP6 P 

Klx{P + Pb,P + Pa+ Pb)D~^{P + Pa + Pb)K^{p + Pa + Pb, P)-D"^(p)] 

= - / tr(75[9,L>(p)][9^L>-^(p)][a,L>(p)][9,L>-^(p)]) (125) 

where the identity 

Dd^D^^ + {d^D)D-'^ = (126) 
and Eqs. (120), (121) and (123) have been used. 

m,n,k,l,j,s,t ^p^y 

P1P2P4 

tr[75L'(pi)L'-^(p2)i^.(p2,P4)i?-'(p4)i^M(p4,Pi)i?-'(pi)] 
= / 5(p„)(5(p,)9fa(^)/tr[75L'(p + p„)L>-i(p + p„+p,) 

K,{p + Pa + P6, p)D-\p)Ki,{p, p + Pa)D-\p + Pa)] 

= / tT(^,[dxDip)][d,D-\p)][d.D{p)][d^D-\p)] 
p 

-^,[dMp)][d.D-\p)][d^D{p)][dxD-\p)]) 

^-2 J tr{^,[d^D-\p)][d.Dip)][dxD-\p)][dMp)]) (127) 
p 

Therefore, summing /„, 1^ and in (122), (124), (125) and (127), and 
using Eqs. (76), (77) and (108), we obtain 



/2 = -4 / tr(^,[d,D-\p)][dMp)][dxD-\p)][d„D 



p 
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(128) 
(129) 



where the operation in (129) produces (128), plus three terms which are 
symmetric in yuz/, ^uA, and /io", respectively, hence neither of these three terms 
contributes to l2- Now we perform the momentum integral in (129) by, first 
removing an infinitesimal ball with center at the origin p = and radius 
e from the Brillouin zone, then evaluating the integral, and finally taking the 
limit e to zero, i.e.. 



7 — 7T lim 

(27r)4 e^O 



df^p dA tr 



<^<\PiJ.\<'^ 



(130) 



Then according to the Gauss theorem, the volume integral over the Brillouin 
zone ( a four dimensional torus due to the periodic boundary conditions ) with 
the ball removed can be expressed as a surface integral on the surface S'g 
of the ball B^, provided that C^{p) is nonzero for e < \p^\ < tt ( i.e., free of 
species doubling ) such that the integrand in (130) is well defined. So, (130) 
becomes 



^—7 lim / 



d s n,,tr 



75 (^2) ^.[^^^2^2 



C2 ^ I 1 + r2C2 



(131) 



where is the /i-th component of the outward normal vector on the surface 
Se- Since we have assumed that C^(p) ^ as p ^ 0, we can set C^(p) = Pfj, 
on the surface S. and obtain 



7 lim / d^s n^tr 

Se 




1 + r'^p'^ 




1 + r^p^ 



(132) 
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where we have used the property 

j £s n^pyfip^) = 5^,^ j d-^s n^p^fip^) (133) 



Finally, we have 



J2 = 



2-7T TT 



^^lim-- — — r [ d(f) [ d92sm92 [ d9i(- sin^ Oicos^ 9i) 







= ^e^^Aa (134) 

This completes the task of evaluating I2 ='H[2 ii{^^RD) ] for R = rl, where 
(134) is one of the main results of this paper. Although (134) has been derived 
for the infinite lattice ( a very large lattice with periodic boundary conditions 
), it is reasonable to expect that it also holds for finite lattices with periodic 
boundary conditions and with even number of sites in each direction, since the 
integral in (130) is essentially a topological invariant quantity. Implications of 
Eq. (134) have been discussed in section 2. 



4 Conclusions and Discussions 



In this paper, we have evaluated the axial anomaly for the Ginsparg- Wilson 
fermion operator D — Df.{l. + rDc)"^. For any chirally symmetric D^. which 
in the free fermion limit is free of species doubling and behaves like i'j^Pfj, as 
p — s> 0, the axial anomaly for U{1) lattice gauge theory with single fermion 
flavor is 

r tr[75D(x, x)] = ■^^eu,^xaFij,,,(x)Fxa{x + fi + u) 

+ higher order and/or non-perturbative terms (135) 

where the field tensor F^y{x) on the lattice is defined in Eq. (10). The FF term 
is r-invariant and has the correct continuum limit. As shown in section 2, the 
higher order and/or non-perturbative terms might have significant impacts 
to the index of D. For smooth background gauge configurations [ e.g., Eq. 
(58)-(61) ] with integer topological charge (62), if the axial anomaly satisfying 
Eq. (50), then the sum of these higher order and/or non-perturbative terms 
over all sites of the lattice is shown to be an integer multiple of the FF term 
and this leads to the emergence of an integer functional, c[D\, which is called 
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topological characteristics of D in ref. [7,8]. In general, c[D] incorporates all 
kinds of contributions from all fermionic modes. Due to the {S, T)-invariance 
of the zero modes [7], we have [ Eq. (64) in section 2 ] 



N_-N+^ c[D] ^ ^f^uXaF^uix)F^,{x + fi + u) (136) 

where c[D] — c[Dc] is invariant for any S and T in the general GW relation, 
Eq. (11). Then for topologically proper D{Dc), i.e., c[D] = 1, the Atiyah- 
Singer index theorem can be realized on a finite ( infinite ) lattice for smooth 
background gauge fields. Now it is obvious that the FF term in Eq. (135) 
must be also (5*, T) -invariant, otherwise it would be contrary to Eq. (136). 
Hence, we have 



1 

- tr{75[(5' + T)D]^^^} = :^^€u,,,xaF,j,i.{x)Fxa{x + pL + u) 

+ higher order and/or nonperturbative terms. 

(137) 

Then one can deduce the following result in the continuum limit. 



^ tT{j,[iS + T)DlJ^c[D] ^e^,xaF^,{x)F^,{x). (138) 

Now the limit ( S, T ^ ) in Eq. (138) can be safely taken. Since the limit 
S*, T — s> is the chiral limit where D = D^. and the GW chiral symmetry 
breaking [ the RHS of Eq. (1) ] is completely turned off, we conclude that the 
GW relation indeed does not play the crucial role to fix the axial anomaly 
of D in the continuum limit. This is in agreement with the conclusion of ref. 
[7]. The crucial point for D to have the correct axial anomaly and to realize 
the Atiyah-Singer index theorem in the continuum limit is the existence of a 
topologically proper which also satisfies the properties mentioned above, 
or in general, the constraints (a)-(e) given in ref. [7]. Then any GW fermion 
operator D constructed by the general solution D = Dc{'S. + RDc)~^ will have 
the desired topological properties. The role of the chiral symmetry breaking 
transformation (5) is to bypass the Nielson-Ninomiya no-go theorem such that 
D can be constructed to be local, free of species doubling and well defined for 
any gauge configurations, while the essential chiral physics of Dc is preserved 
under this transformation. Therefore, for practical computations on a finite 
lattice ( with finite lattice spacings ), one must keep {S,T) finite as well as 
using a topologically proper D{Dc) such that the axial anomaly could agree 
with the Chern-Pontryagin density in continuum, though the index is equal 
to the topologcial charge for any [S, T). 
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A 



In this appendix, we explicitly show that 

'H'[eai3^sFa(3ix)F^six + a + p)]= 'H'[eai3^sFai3ix)F^s{x)] = Sei^^xa (A.l) 
where H' is defined as 



n'[0]=J2x,yx 



x,y 



[0{x)] 



(A.2) 



which is similar to H, defined in Eq. (47) but without imposing the gauge fields 
to zero after the differentiations with respect to the gauge fields. 

First, wc derive the field tensor for the U{1) gauge theory on the lattice. A 
plaquette on the fi — i> plane is defined as 



P^.u{x) = U^{x)U,{x + ix)Ul{x + v)Ul{x) 

= exp|iea [A^{x + jl) - Ai,(x) - + u) + 



and its expansion up to is 



1 + iea 
1 

'2 



2 2 

— e a 



A^{x + A^{x) - A^{x + i>) + A^{x 
A^{x + 11)- A^{x) - A^{x + z>) + A^{x) 



~ 1 + iea 



dfiAn diiA^ 



-^2 4 

e a 

2 



d^A^-d^A^' + 0{e\a^^ 



Then the real part of the sum of all plaquettes, i.e.. 



^EERe[i-p,.(x)], 



goes to 



(A.3) 



(A.4) 



(A.5) 



a^ \di,A^ - dyA^ 



(A.6) 



in the continuum limit, thus agrees with the action of QED. Hence, we can 
identify the field tensor on the lattice to be 
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F^,{x) = - [A,{x + fi)- A,{x) - A^{x + z>) + A^{x)] (A.7) 

(X 

= ^log[P^.(a;)] (A.8) 

We note that on a finite lattice with periodic boundary conditions, and for 
background gauge fields with nonzero topological charge, some of the link 
variables at the boundary need modifications [ Eqs. (59) and (61) ] such that 
the field tensors are continuous on the torus. 

Using Eqs. (A. 2) and (A.7), we obtain 



T-C' [eap^sFap {x) F^s {x)] 

^Afs{x - a)As{x - 7) - Aa{x - i3)As{x - 7) 
-Ap{x - a)A^{x -S)+Aa{x- P)A^{x - S) ^ 

i^auX/j. ~ ^uaXn ~ ^cvfxX + ^v<Ttx\)^v{x — (t)A^(x — A) 
+ (fA/io-i/ ~ ^nXav ~ ^Xixvcr + ^nXva)-^iJ,{x — X)Ai,{x — (t) | 

= 8e^,Aa (A.9) 

where in the first and the second equalities, only those terms which have 
nonzero contributions are retained. 

For Liischer's topological charge density having the second field tensor located 
at the site x + ji + i) rather than at we obtain 



[ Ap{x + a)As{x + a + /9 + 7) - Ap{x)As{x + a + /9 + 7) 

-Aa{x + (5)As{x + a + /3 + 7) + Ao,{x)As{x + a + + 7) 
-Ai3{x + a)A^{x + a + /9 + 5) + Ai3{x)A^{x + a + ^9 + 5) 
+A^{x + ^)A^{x + a + ;9 + 5) - A^{x)A^{x + a + ^ + 5) 

^S'^'^^^M.(0)M,(y){ 

^ai'XtiAy{x + a)A^,{x + (T + i> + A) + ea^xuA^,{x + a)Ay{x + a + ji + X) 
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- ^aiJ,XuAi^{x)A,,{x + a + fi + \) - exij,auAij,{x)Ai,{x + X + fi + a) 

- ^uaXix^uix + ^)Ai^{x + i> + a + X) - e^^xvA^{x + a)A^{x + /i + i> + A) 
+ ^iJXT\uA^{x)A„{x + /t + a- + A) + efj,xa^A^{x)A„{x + /t + A + ct) 

- ^av^l\Au{x + a)A^{x + o" + i> + A) - e^^j,uxA^{x + a)A,,{x + + /t + A) 
+ (■\,,uaA^{x)A^{x + A + /t + 0-) + e^f,yxA^{x)A^{x + a- + /i + A) 

+ e,.aMA^!/(a; + a)Afj_{x + i> + a + A) + e^^^^A^^l^ + + fi + a + \) 

- e^,xu^A^{x)Ay{x + /i + A + a) - e^ayxA^{x)Ay{x + /i + a + A)| 




Among the eight hnes of expressions at the second equaUty, the second hne 
vanishes due the cancellation of its two terms, and the same happens to the 
fourth, the sixth and the eighth lines. Then the remaining four lines add up 
to yield the final result. This completes the proof of the identity (A.l). 



In this appendix we derive some useful properties of the kernel of the 
vector current 



These properties [ Eqs. (B.IO) - (B.12) ] are given in the appendix of Ginsparg 
and Wilson's original paper [1]. Here we present our derivation in details and 
correct a minor misprint in ref. [1]. As usual, the divergence of the vector 
current is extracted from the change of the action under an infinitesimal local 
transformation at the site n, 

lK^1pn + On'4'n 

with the prescription 



B 




m,l 



then we obtain 



(B.2) 



m 
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where 9^J|u(n) is defined by the backward difference 

d^J,{n) = Y: [JM - Mn - A) ]. (B.3) 

Since we will turn off the gauge field after performing the differentiations in 
Eq. (73), we only need to derive all properties of in the free field hmit 
where the action is translational invariant Dj^^ — -C^m-n- Then Eq. (B.2) can 
be rewritten as 

d^Jnin) = Y^{'4)n+lDl-4^n - i>nDli>n-l)- (B.4) 

I 

and Eq. (B.l) as 

•^m(^) = Yl '^n+mK^{m, m - l)'ll;n+m-l- (B.5) 

m,l 

To construct Kfj_{m,m — I) such that Eq. (B.4) can be reproduced with Eqs. 
(B.3) and (B.5), the authors of ref . [1] set 

Ku,(m, m-l) ^ Lfj,{m, I) Di = sign(y /^(m, /) A- (B.6) 

where /^(m, /) is equal to the fraction of the shortest length paths from to I 
which pass through the hnk from m — fi to m. It is straightforward to verify 
that this definition of leads to Eq. (B.4). Using Eq. (B.5) and Eq. (B.3), 
we obtain 

li,m,,l 

-iJn-fX+mK^{m, m - l)lpn-iX+m-l] (B.7) 

The number of shortest length paths from n + m — I to n + m is 





1 + 1 




+ 




+ 1 


h\)\ 




^ll 


I 


Ik 


I 




1^4 


I 



For a given set of positive intgers (81,82,83,84), and for all I with — 
Su,i^ = I, • • • ,4:, then for > 0, each one of these paths passing through the 
hnk from n + m — fi to n + m contributes ■^^■xpn+m.Diipn+m-i to J^{n); while for 
^/i ^ 0, a path passing through the link from n + m to n + m — fi contributes 
— ^^ipn+mDiipn+m-i to J^{n). Hcucc for cach /, the contribution of each link 
to the RHS of Eq. (B.7) is 
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'irri'^Pn+mDltpn+m-l ~ '^n-fi+mDiXpn-jl+m-l) ■ 

Adding their contributions along each shortest length path would cancel in 
pairs except the boundary terms 



^{tpn+lDltpn - 1pnDl1pn-l)- 

The sum over all shortest length paths then cancels the factor Thus Eq. 
(B.4) is reproduced. 

Next we prove two identities which are essential for deriving Eqs. (B.IO) - 
(B.12), 



Y,L,{m,l)=l,. (B.8) 

m 

Y,m^L,{m,l)J-^^^^^^±^. (B.9) 

m ^ 

The summation over sites m is defined as 



Iv la l-X 

E-E E E E 

m m/j, m,^=0 ma=0 mx=0 

where the upper and lower limits of depend on the sign of l^. For > 0, 
the summation of is from 1 to l/^, while for 1^ < 0, from to + 1. 



E- 



E'^^^i if > 
T.'^-^Lo if < 



To prove the first identity, we observe that ( see Fig. 1 ) all shortest length 
paths from to I must go through one of the links pointing in the jl di- 
rection with a fixed m^, i.e., they are all perpendicular to the hyperplane 
with fixed m^^. Therefore holding fixed and summing over all other in- 
dices ( my^iTifj.mx ) of the fraction f^{m^l) [ defined in Eq. (B.6) ] is equal 
to summing all probabilities for all shortest length paths going through the 
hyperplane and hence it must equal to one. Since there are hyperplanes 
between and Z^, after summing over m^, we obtain Eq. (B.8), 
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m- \i 












m 







1^ 



Fig. B.l. The portion of the lattice containing the shortest length paths between 

and I. Only two directions ( fi and u ) arc shown on the plane while the other 
two directions ( a and A ) are orthogonal to the plane. The solid line with arrows 
denotes the projection of one of the shortest length paths onto the plane. 

Iv la l\ 

J2 J2 sign(y/^(m,0 = ^sign(y = 
For the second identity, we first prove the case — ly, 



ll/ 1(7 l\ 

HY. I] rn^^/xK/) = l^m/, sign(/^) 

[1 + 2 + --- + Q if > 

[-1-2 + if l^.<0 

2 

For z/ 7^ /i, by symmetry, the fraction of shortest length paths going through 
those links with fixed and my is ( see Fig. 1 ), 



la h 

Y Y 

mcr=0 m\=0 



1 



|/J + 1 



Then 



li/ la l\ lu 'Y 

YY Y Y myL^{m,l) = YY sign(/^) 

mu=0 mcy=0 m\=0 "ip mi,=0 

This completes the proof of the second identity. 



With these two identities, we proceed to derive some useful properties of 
in momentum space 



l,m 
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l,m 

Since we only need the derivatives of with respect to p'^^ evaluated at p' — 0, 
we can expand the above equation up to first order of p' 



K^{p,P + p') « E + ^P'{m - l)]L^{m, l)Di 

l,m 



where D{p) = J2i^~^^^Di and Eqs. (B.8), (B.9) have been used. Therefore, we 
have 



K^{p,p)=td^D{p) (B.IO) 

and 

dlK^{p,P + P%'=o = {^d^d, - ^S,,d^^ D{p) (B.ll) 
Similarly, we obtain 



K^^ip + p\p) = E e-^(^+^')'"e*^('"-')i^^(m, m - I) 

l,m 
l,m 

~ X] - ip'm)Kn{m, m - I) 

l,m 



2 

=«.(i + ^-f)^(p) 



and this leads to 



KK^{p + p',p)\p'=o = {^d,d, + ^S,,d,^ D{p) (B.12) 

We note that the operator in the second term of Eq. (B.ll) is missed in 
Eqs. (37) and (A6) of Ginsparg and Wilson's original paper [1]. Equations 
(B.10)-(B.12) are used in our derivation of chiral anomaly in section 3. 
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c 



In this appendix we prove the second equahty of Eq. (92) [ or Eq. (93) ] which 
has played an important role in our derivation of axial anomaly in section 3. 
The vector current in Eq. (92) is 

Jfi{k,U) = J2^k+iK^{k,i,j;U)il:k+j = i^^m j~. ,,s Dmn{U)lpn (C.l) 

i,j m,n (^^iJ,yi^) 



An explicit realization of the kernel K^{k,i,j; U) is given in ref. [1] as 

K^{k, i, j- U) = sign((z - j)^)f^{k + i,k + j)Dk+i,k+j{U) (C.2) 

where /^(/c + i,k + j) is the fraction of the shortest length paths from k + j 
to k + i which pass through the link (A:, k + fi). Note that the vector current 
J^(A;, U) satisfying Eq. (36) is not unique. But the kernel defined in Eq. (C.2) 
leads to the second equahty of Eq. (C.l). First we note that the ordered 
product of link variables along one of the shortest length paths from n to m, 
say, P, is 

Up{m,n) =]lU{m,s)---U{t,n) (C.3) 
p 

where U {m, s) denotes the link variable pointing from m to s with the usual 
convention 

U{m,m + fi) = f/^(m) = exp[meA^(m)]. 

Then Up{rn,n) enters the action A = J2m,n'4^mDmn{U)ipn via the following 
gauge invariant product 

'ipm^{m,n)Up{m,n)ijjn (C.4) 

where r(m, n) is a matrix in the Dirac space and its explicit form is irrelevant 
to our present discussion. The normalized sum of (C.4) over all shortest length 
paths from n to m is equal to the term '0j„L>j„„(C/)'0„ in the action, i.e., 

D,^n{U) = r{m,n)^YlUp{m,n) (C.5) 

where Ni is the total number of shortest length paths from n to m, 

(I/1I + I/2I + I/3I + I/4I)! , 
= 1, II I, II I, I, I, I, > l^m-n (C.6) 
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Then using the following simple identity 

6 



U^{n) = i5i^JknUf^in) (C.7) 



it is straightforward to obtain that, if Up{m, n) contains the link U^{k) ( or its 
hermitian conjugate ), then the derivative of DmniU) with respect to A^{k) 
yields i Up{m,n) times a sign factor which depends on the relative positions 
of m and n ( i.e. —1 for > but +1 for < [ see Eq. (C.4) ] ); and 
the derivative is zero if Up{m,n) does not contain the link U^{k). That is 



-i sign{{m - n)i^)Up{m,n), if C/^(/c) e Up{m,n) 
0, otherwise 



Now multiplying both sides of Eq. (C.8) by r(m, n), summing over all shortest 
length paths from n to m, and dividing by Ni, we obtain 



S 1 
(fe) ^"""*^^^ " sign((m - n)^)r(m, ^)]^ ("^' ^) 

= sign((m - n)i,)fi^{m,n)r{m,n)^^Up'{m,n) (C.9) 

where Eq. (C.5) has been used on the LHS. The summation P' on the RHS 
denotes the sum over the shortest length paths between m and n which pass 
through the link {k, k + fi), and the total number of these paths is denoted by 
Ni, and f^{m,n) = NjNi is the fraction of the shortest length paths from n 
to m which pass through the link (/c, k+fi). Now we sandwich both sides of Eq. 
(C.9) by %l)m and ipn and sum over m and n. Then on the RHS of the resulting 
equation, we can write m — k + i and n — k + j and replace summations 
over m and n by summations over i and j, since only those shortest length 
paths from n to m which pass through the link [k, k + p.) can have nonzero 
contributuion, finally we have 



^ ^m-Fl-7TT^mn(^)^r, 



^J2'^k+i sign((i - j%) fu,{k + i, k + j) Dk+i,k+j{U) i/^k+j (C.IO) 
^Y.i^k+iKf,{k,i,j;U)i/:k+j (C.ll) 
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where 



Dk+i,k+j{U) = T{k + t,k + Up'ik + t,k+j). (C.12) 



This completes the proof of the second equahty of Eq. (92). 
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